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The classical Volterra predator-prey model concerns the first order planar differ- 
ential system 



where a, b, c, d > are constant coefficients. The study of system (Eq) is confined 
to the open first quadrant (Mq )^ of the plane, since x(t) > and y{t) > represent 
the size (number of individuals or density) of the prey and the predator populations, 
respectively. Such model was proposed by Vito Volterra in 1926 answering D'Ancona's 
question about the percentage of selachians and food fish caught in the northern 
Adriatic Sea during a range of years covering the period of the World War I (see 
[H for a more detailed historical account). 

System (Eq) is conservative and its phase-portrait is that of a global center at the 
point 



surrounded by periodic orbits (run in the counterclockwise sense) which are the level 
lines of the first integral 



1. Introduction and main results 





(1.1) 



So{x, y) := dx — clogx + by — a\ogy, 



(1.2) 
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that we'll call "energy" in analogy to mechanical systems. The choice of the sign 
in the definition of the first integral implies that SQ{x,y) achieves a strict absolute 
minimum at the point Pq . 

According to Volterra's analysis of (Eq), the average of a periodic solution {x(t),y{t)), 
evaluated over a time-interval corresponding to its natural period, coincides with the 
coordinates of the point Pq • 

In order to include the effects of fishing in the model, one can suppose that, during 
the harvesting time, both the prey and the predator populations are reduced at a 
rate proportional to the size of the population itself. This assumption leads to the 
new system 

/ x' = x{a^-by) , . 

\ y' = y{-c, + dx), ^"^^^ 

where 

:= a — jj, and : = c + /i 

are the modified growth coefficients which take into account the fishing rates —ixx{t) 
and —fiy{t), respectively. The parameter fi is assumed to be positive but small enough 
{fi < a) in order to prevent the extinction of the populations. System {Efj_) has the 
same form like (Eq) and therefore its phase-portrait is that of a global center at 

The periodic orbits surrounding are the level lines of the first integral 

Sf,{x,y) := dx - c^^logx + by - a^logy. (1.4) 

The coordinates of P^ coincide with the average values of the prey and the predator 
populations under the effect of fishing (see Figure 1). A comparison between the coor- 
dinates of Po and P^ motivates the conclusion {Volterra's principle) that a moderate 
harvesting has a favorable effect for the prey population [1]. 

The Volterra's system (also called Lotka- Volterra's with reference to the work of 
Alfred Lotka, who first in 1920 used the same system as a description for certain 
chemical reactions) has been sometimes criticized by ecologists and biologists, who 
refused to consider such a model as accurate (see [3] for a discussion about this topic), 
and, in the course of the years, many variants of it have been proposed. Volterra's 
himself modified system (Eq) in [HI B5] , by replacing the Malthusian growth rates 
with logistic terms of Verhulst type (see also [26]). In order to incorporate the effects 
of a cyclic environment, periodic coefficients have been introduced both in the basic 
model and in its variants. The past forty years have witnessed a growing interest 
in such kind of models and several results have been obtained about the existence, 
multiplicity and stability of periodic solutions for Lotka- Volterra type predator-prey 
models with periodic coefficients [Il[2l[3l[6l[7l[g[IIl[l2l[l3l[T5l[l6l[ni[23l[2l 
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Figure 1 . in this picture we show some periodic orbits of the Volterra system 
(Eq) with center at P = Pq as well as of the perturbed system (i?^) with center at 
Q = Pf^ (for a certain /i g ]0, a[ ). 

If we take the Volterra's original model (Eo) and assume a seasonal effect on the 
coefficients, we are led to consider a new system of the form 



where a(-), b{-), c(-), d{-) : M ^ M are periodic functions with a common period T > 0. 
In such a framework, it is natural to look for harmonic (i.e. T-periodic) or m-th 
order subharmonic solutions (i.e. mT-periodic, for some integer m > 2, with mT the 
minimal period in the set {jT : j = 1, 2, . . . }) having range in the open first quadrant 
[positive solutions). With this respect, we have the following theorem which can be 
derived as a corollary of some results in [T3] dealing with certain classes of time- 
periodic Kolmogorov systems. In Theorem 11.11 below, as well as in the other results 
of this paper, solutions are meant in the Caratheodory sense, that is, (x(t),y(t)) is 
absolutely continuous and satisfies system (E) for almost every t G M. Of course, such 
solutions are of class if the coefficients are continuous. 

Theorem 1.1. Suppose that b{-) and d{-) are continuous functions such that 

h{t) > 0, d{t) > 0, yte [0,T] 
and let a(-),c(-) G L^{[0,T]) be such that 



x' = x{a{t) — b{t)y) 
y' = y{-c{t) + d{t)x) 



a := — 




a{t) dt > 0, c 




c(t) dt > 0. 



4 



M. PIREDDU, F. ZANOLIN 



Then the following conclusions hold: 

(ci) System (E) has at least one positive T -periodic solution; 

(62) There exists an index m* > 2 such that, for every m > m*, there are at 
least two suhharmonic solutions {xi^rn,yi,m) and {x2,m,y2,m) of order m for 
(E) which do not belong to the same periodicity class and satisfy 

lim^^oo (mill Xj,™) = lim^^oo (min iji^rn) =0 (i = I 2) 

lim^^oo {max Xi^rn) = lim^^oo (maxyj^m) = +00 ^ ' 



For a proof and other details, see [13] . We refer to [12] for detailed information about 
the subharmonic solutions of system (E) and to [H, [231 123] for results about the stabil- 
ity and the number of solutions. See also [H [HI [13] for more general conditions on the 
coefficients ensuring a priori bounds and existence of T-periodic positive solutions. 

Let us come back for a moment to the original Volterra system with constant 
coefficients and suppose that the interaction between the two populations is governed 
by system {Eq) for a certain period of the season (corresponding to a time-interval 
of length To) and by system [E^) for the rest of the time (corresponding to a time- 
interval of length r^). Assume also that such alternation between (Eq) and (E^j) 
occurs in a periodic fashion, so that 

T := ro + r^ 

is the period of the season. In other terms, we consider at first system (Eq) for 
t e [0, ro[. Next we switch to system {E^) at time tq and assume that {Efj_) rules the 
dynamics for t G [ro,T[. Finally, we suppose that we switch again to system (Eq) at 
time t = T and repeat the cycle with T-periodicity. 

Such two-state alternating behavior can be equivalently described in terms of equa- 
tion {E), by assuming 

a(t) = af,(t) 
c{t) = c,it) 

as well as 

b{t) = b, d{t) = d, 

with a, b, c, d positive constants and /i a parameter with < n < a. Hence we can 
now consider system 

j x' = x{a^{t)-hy) , . 

\ y' = y(-c^(t)+dx), > 

where the piecewise constant functions and are supposed to be extended to 
the whole real line by T-periodicity. Clearly, Theorem 11.11 holds for equation [E*), 
ensuring the existence of at least one positive T-periodic solution and m-th order 
subharmonics of any sufficiently large order. 



(a for < t < ro 

a — n for ro < t < T , 

J c for < t < ro 

c + /i for ro < t < T , 
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It is our aim now to prove that system (E*) generates far richer dynamics. Indeed, 
we'll show the presence of a topological horseshoe for the Poincare map 

: (M+)2 ^ (M+)', <Piz):=C{T,z), 

where C(')^) = ^)) is the solution of (E*) starting from z = {xo,yo) G 

(M(]')^ at the time t = 0. As a consequence, all the complexity which is associated to 
an horseshoe geometry (like, for instance, a semiconjugation to the Bernoulli shift, 
sensitivity with respect to initial conditions, positive topological entropy, a compact 
set containing a dense subset of periodic points) will be guaranteed. To this goal, we 
apply recent developments [32] which connect the analysis of certain planar ODEs 
to the theory of linked twist maps. With such a term, one usually designates some 
geometric configurations characterized by the alternation of two planar homeomor- 
phisms (or diffeomorphisms) which twist two overlapping annuli. More precisely, we 
have two annular regions A and B which cross in two disjoint topological rectangles 
TZi and 7^2 • Each annulus is turned onto itself by a homeomorphism which leaves 
the boundaries of the annulus invariant. Both the maps act in their domain so that 
a twist effect is produced. This happens, for instance, when the angular speed is 
monotone with respect to the radius. Under certain assumptions, it is possible to 
prove the existence of a Smale horseshoe inside TZi {i = 1,2) [TU]. Linked twist maps 
(LTMs) were originally studied in the 80s by Devaney [10], Burton and Easton [S] and 
Przytycki [3SIEZ!- As observed in [TD], such maps naturally appear in mathematical 
models for particle motions in a magnetic field and in differential geometry. Geomet- 
rical configurations related to LTMs can be also found in the restricted three-body 
problem [271 PP- 90-94]. In more recent years significant applications of LTMs have 
been performed in the area of fluid mixing (see, for instance, [l2l HHl H9]). 

With the aid of Figure 2 we try now to explain how to show the presence of a 
horseshoe-type geometry for the switching system (E*). As a first step, we take two 
closed overlapping annuli made up by level lines of the first integrals associated to 
system [Eq) and (Ef^), respectively. In particular, the inner and outer boundaries 
of each annulus are closed trajectories surrounding the equilibrium point (P = Pq 
for system (Eq) and Q = for system {E^)). Such annuli (that we call from now 
on Ap and Aq) intersect in two compact disjoint rectangular sets TZi and 7^2 • The 
order in which we decide to name the two regions TZi and 7^2 is completely arbitrary. 
Whenever we enter in a setting like that described in Figure 2, we say that the annuli 
Ap and Aq are linked together. Technical conditions on the energy level lines defining 
Ap and Aq , which ensure the linking condition, are presented in Section [2l 

As a second step, we give an "orientation" to TZi (for i = 1,2) by selecting, in the 
boundary, two disjoint compact arcs (i.e. homeomorphic images of the compact unit 
interval of the real hne) that we denote by 'R-~ i^^^ and Ti'^ j-ig^t and call the left and 
the right sides of dTZi . We also set 

■~ '^i,left ^ '^i, right ■ 

The closure of dTZi \ TZ~ is denoted by TZf . It consists of two disjoint components 
which are compact arcs that we name the down and up sides of dTZi (for a precise 
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Figure 2. The two annular regions Ap and Aq (centered at P and Q, respec- 
tively) are linked together. We have drawn with a darker color the two rectangular 
sets TZi and TZ2 where they meet. 

definition of oriented rectangle, see Section [3]). In tlie specific example of Figure 2, 
we orientate TZi and 7^2 as follows. We take as 71^ the intersection of TZi with the 
inner and outer boundaries of and as 7^2 intersection of 7^2 with the inner 
and outer boundaries of Aq . The way in which we choose to name (as left/right) the 
two components of 7l~ is inessential for the rest of the discussion. Just to fix the 
ideas, let's say that we call "left" the component of 7^]" which is closer to P and the 
component of 7^^ which is closer to the equilibrium point Q (of course, the "right" 
components will be the other ones). 

As a third step, we observe that the Poincare map associated to (E*) can be 
decomposed as 

= 0^ O 00 , 

where 0o is the Poincare map of system {Eq) on the time-interval [0, tq] and 0^ is the 
Poincare map for (E'^) on the time-interval [0,r^] = [0,T — tq]. 
Consider also a path 7 : [0, 1] TZi with 7(0) G 7^|f;e/t ^ right ■ -^^ we'll 

see in Section |2l the points of 71^ i^j-^ move faster than those belonging to 71]^ ^^^^^^ 
under the action of system {Eq). Hence, for a choice of the first switching time tq 
large enough, it is possible to make the path 

[0,1] 9^^^0o(7W) 
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turn in a spiral-like fashion inside the annulus Ap and cross at least twice the rect- 
angular region 7^2 from 7^^ ^^^^ to 7^^ . Thus we can select two sub-intervals of 
[0, 1] such that 0o ° 7 restricted to each of these intervals is a path contained in 7^2 
which connects the two components of 7^.2 • 

Now, we observe that the points of 7^2 left ^ove faster than those belonging to ^^gj^^ 
under the action of system (-E^). Therefore, we can repeat the same argument as above 
and conclude that for a suitable choice of = T — tq large enough we can transform, 
via , any path in 7^2 joining the two components of 7^^ onto a path which crosses 
at least once TZi from 71^ ^^^^ to 1Z~[ ^^^^^^ . 

As a final step, we complete our proof about the existence of chaotic-like dynamics 
by applying a topological lemma that we recall in Section [3] as Lemma 13.11 for the 
reader's convenience. 

In conclusion, our main result can be stated as follows. 

Theorem 1.2. For any choice of positive constants a,b,c,d,fi with fi < a and for 

every pair {Ap,Aq) of annuli linked together, the following conclusion holds: 

For every integer m > 2 there exist two positive constants a and /5, such that for each 

ro > a and > /?, 

the Poincare map associated to system {E*) induces chaotic dynamics on m symbols 
in TZi and in 7^2 • 

In view of our result, one could conclude that complex dynamics were already 
hidden in Volterra's work ^4J, since linked twist maps (on long periods) appear as a 
consequence of the monotonicity of the period map and of two Volterra's principles: 

• "Se si cerca di distruggere uniformemente e proporzionalmente al loro numero 
gli individui delle due specie, cresce la media del numero di individui della 
specie mangiata e diminuisce quella degli individui della specie mangiante" ; 
(If one tries to destroy uniformly and proportionally to their number the individuals of the 
two species, the average of the number of individuals of the eaten species increases and the 
one of the eating species decreases); 

• "Se si distruggono contemporaneamente e uniformemente individui delle due 
specie, cresce il rapporto delVampiezza della fluttuazione della specie mangiata 
all'ampiezza della fluttuazione della specie mangiante". 

(If one destroys at the same time and uniformly the individuals of the two species, the 
ratio between the amplitude of the fluctuation of the eaten species and the amplitude of 
the fluctuation of the eating species increases). 

Indeed, the first principle says that the position of the center, around which the annuli 
can be constructed, varies according to the strength of the fishing, while the second 
principle implies that the shapes of the annuli are suitable for a linking (see [44, fig. 
6]). The twist conditions on the boundaries for the associated Poincare maps come 
from an analysis on the periods of the orbits (which was carried on by Volterra in the 
limit case, i.e. for orbits near the equilibrium points). 
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In order to clarify the meaning of Theorem 11.21 we introduce the precise concept 
of chaotic dynamics that we consider in this work. Our definition is a modification 
of the corresponding one in [20] and abstracts the usual interpretation of chaos as 
the possibility of realizing any coin-flipping sequence, by giving also a special em- 
phasis to the presence of periodic orbits. Definitions presenting similar features have 
been considered by several authors dealing with nonautonomous ODEs with periodic 
coefficients [H HH [50], as well as in abstract theorems about periodic points and 
chaotic-like dynamics in metric spaces ^0] . We refer to Section [3| for a more detailed 
discussion about the kind of complex dynamics involved in Definition 11.11 below. 

Definition 1.1. Let Z he a metric space and ip : Z D —>■ Z he a map. Let also 
1^ ^ he a nonempty set and m > 2 be an integer. We say that ip induces chaotic 
dynamics on m symbols in D if there exist m pairwise disjoint (nonempty) compact 
sets 

/Ci , /C2 , . . . , /Cm CD, 
such that, for each two-sided sequence of symbols 

(Si)i6z e Em := {1, . . . ,m}^ , 
there exists a two-sided sequence of points 

with 

Wi G Kg. and t^j+i = 4'(wi), Vi G Z. (1.5) 

Moreover, if (sj)igz is a /c-periodic sequence (that is, Si+k = Si,Wi G Z) for some 
k > 1, then there exists a fc-periodic sequence {wi)i^z satisfying (11.51) . When we 
wish to put in evidence the role of the sets /C^'s, we'll also say that induces chaotic 
dynamics on m symbols in the set T> relatively to (/Ci, . . . , /Cm)- 

From this definition it follows that if ip is continuous and injective on the set 

m 

/C := IJ/Cj CD, 

then there exists a nonempty compact set 

JC /C, 

for which the following properties are fulfilled: 
• X is invariant for ijj (i.e. ip{T) = X); 



Such assumptions are fulfilled in our application: indeed, in Thcorcm ll.21 ^ — (/} is a, homeomor- 
phism, being the Poincare map associated to (E*). 
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is semi conjugate to the Bernoulli shift on m symbols, that is, there exists 
a continuous map g oil onto such that the diagram 



X 



y K V 



;i.6) 



commutes, where a is the Bernoulli shift on m symbols ( 
the homeomorphism defined by a{{si)i) := Vi G Z); 

• The counterimage g~^{s) C X of every /c-periodic sequence s = (sj)jgz G 
contains at least one /c-periodic point. 

For a proof, see Lemma 13.21 in Section [3l From the above properties it also follows 
that 

^top(^lx) > htopia) = log(m), 
where htop is the topological entropy [17] . Moreover, according to [I8| Lemma 4], 
there exists a (nonempty) compact invariant set X' C X such that has sensitive 
dependence on initial conditions, i.e. 3(5 > : Vw G X' there is a sequence Wi of 
points in X' such that Wi ^ w and for each i G N there exists m = m(i) with 
dist(V^"^(wi),^™(w)) > ^. 

As remarked in [32j, if we look at Definition II. II and its consequences in the context 
of concrete examples of ODEs (for instance when ip turns out to be the Poincare map), 
condition (11.51) may be sometimes interpreted in terms of the oscillatory behavior of 
the solutions. Such situation occurred in [211 [SS] and takes place also for system 
{E*). Indeed, from the proof of Theorem 11.21 one sees that it is possible to provide 
more precise conclusions in the statement of our main result. Namely, the following 
additional properties can be proved: 
For every splitting of the integer m > 2 as 

m = mim2, with mi,m2 GN, 

there exist nonnegative integers ni , K2 , with ki = Ki{rQ,mi) and k,2 = '^2(^5 "^2); 
such that, for each two-sided sequence of symbols 

s = (sj)iez = {Pt, qi)i& G {1, . . . , mi}^ x {1, . . . , mg}^ , 
there exists a solution 

Cs(-) = (^s(-),2/s(-)) 

of {E*) with Cs(0) G Til such that Qsif) crosses 7^2 exactly ni + pi times for t G 
]ir, rQ + iT[ and crosses IZi exactly K2 + q-i times for t G ]ro + ^T, (^ + 1)T[ . Moreover, 
if (si)jez = {Pi,(li)i& is a periodic sequence, i.e. Si+k = Si for some k > 1, then 
Ut + kT) = Ut),WteR. 

Note that also the splitting m = mi 777,2 with mi = m and 7772 = 1 (or mi = 1 and 
7772 = "o^) is allowed. For the precise connection between the role of the compact sets 
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A^j's in Definition 11.11 and the oscillatory behavior of the solutions, we refer to the 
proof of Theorem 11.21 in the next section. 

The constants a and (5 which represent the lower bounds for tq and can be esti- 
mated in terms of m\ and and other geometric parameters, like the fundamental 
periods of the orbits bounding the linked annuli (see (12.11) and (12.51) ). 

We end this introductory section with a few observations about our main result. 

First of all, we note that, according to Theorem ll.2t there is an abundance of chaotic 
regimes for system (-E*), provided that the time-interval lengths rg and (and, 
consequently, the period T) are sufficiently large. Indeed, we are able to prove the 
existence of chaotic invariant sets inside each intersection of two annular regions linked 
together. One could conjecture the presence of Smale horseshoes contained in such 
intersections, like in the classical case of the linked twist maps with circular annuli as 
domains |10j. On the other hand, in our approach, which is purely topological (like 
similar ones proposed in [TU [TSl EI]), we just have to check a twist hypothesis on 
the boundary, without the need of verifying any hyperbolicity condition. Hence, our 
technique allows to detect the presence of chaotic features by means of elementary 
tools. This, of course, does not prevent the possibility of a further deeper analysis 
using more complex computations. We also observe that our result is stable with 
respect to small perturbations of the coefficients. Indeed, as it will be clear from the 
proof, whenever tq > a and > jS are chosen in order to achieve the conclusion 
of Theorem II. 2[ it follows that there exists a constant e > such that Theorem 11.21 
applies to equation (E) too, provided that 



Here, the T-periodic coefficients may be in L^([0,T]) or even continuous or smooth 
functions (possibly of class C°°). 

A further remark concerns the fact that, in our model, we have assumed that the 
harvesting period starts and ends for both the species at the same moment. With 
this respect one could face a more general situation where some phase-shift between 
the two harvesting intervals occurs. Such cases have been already explored in some 
biological models, mostly from a numerical point of view. See for an example 
dealing with competing species and [38] for a predator-prey system. If we assume a 
phase-shift in the periodic coefficients, that is, if we consider 



for some < 61,62 < T, and we also suppose that the length tq of the time-intervals 
without harvesting may be different for the two species (say tq = e]0,T[ in 
the definition of and vq = g]0,T[ in the definition of c^), then the geometry 
of our problem turns out to be a combination of linked twist maps on two, three 
or four annuli (which are mutually linked together). In this manner, we increase 




a(t) := d^(t — 61) and c(t) := c^(t — 62) 



CHAOTIC DYNAMICS 11 

the possibility of chaotic configurations, provided that the system is subject to the 
different regimes for sufficiently long time. For a pictorial comment, see Figure 3 
where all the possible links among four annuli are realized. 




Figure 3. We have depicted four linked annular regions bounded by energy level 
lines corresponding to Volterra systems with centers a,t P = Pq, P' = (c/d, 
Q = Pfj, and Q' = {c^/d, a/b), by putting in evidence the regions of mutual inter- 
section, where it is possible to locate the chaotic invariant sets. 



As a final observation, we notice that our approach can be also applied (modulo more 
technicalities) to those time-periodic planar Kolmogorov systems [21j 

x' = X{t,x,y), y' = Y{t,x,y), 

which possess dynamical features similar to the ones described above for the Volterra's. 
Investigations in this direction will be pursued elsewhere. 

The paper is organized as follows. In Section [2] we provide the details for the 
proof of Theorem 11.21 in the sense that, following the argument described above, we 
justify all the steps by means of some technical estimates. In Section [3] we recall the 
topological tools, the corresponding notation and the abstract theorems which are 
employed along the paper. We conclude the article with a brief discussion on chaotic 
dynamics in the sense of Definition 11.11 
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2. Technical estimates and proof of the main result 
Let us consider system (Eq) and let 

i> Xo-= £q{Pq) = ram{£o{x,y) : x > 0,y > 0}. 

The level line 

To{i) :={(x,y)G(M+)' : So{x,y) = i} 
is a closed orbit (surrounding Pq) which is run in the counterclockwise sense, complet- 
ing one turn in a fundamental period that we denote by To{i). According to classical 
results on the period of the Lotka-Volterra system [391 SS], we know that the map 

is strictly increasing with ro(+oo) = +oo and satisfies 

lim To{i) =To :-- 



l^Xo V 

Similarly, if we consider system (E^) with < yU < a, for 

h> Xt,-= Sf^iPfi) = inm{£f,{x,y) : x>0,y>0}, 

we denote by T^{h) the minimal period associated to the orbit 

T^{h) :={(x,i/)G(M+)' : £,ix,y) = h}. 

Also in this case, we have that the map h i— r^(/i) is strictly increasing with r^(+oo) 
+00 and 

lim r^(/i) = 



Before giving the details for the proof of our main result, we describe conditions 
on the energy level lines of two annuli Ap and Aq, centered at P = Pq = f) 
and Q = = (^^, ^-^) , respectively, sufficient to ensure that they are linked 
together. With this respect, we have to consider the intersections among the closed 
orbits around the two equilibria and the straight line r passing through the points P 
and Q, whose equation is by + dx — a — c = 0. We introduce an orientation on such 
line by defining an order among its points. More precisely, we set A ^ P (resp. 
A -< P) if and only if < xb (resp. xa < xb), where A = {xa, yA), B = {xb, ys)- In 
this manner, the order on r is that inherited from the oriented x-axis, by projecting 
the points of r onto the abscissa. Assume now we have two closed orbits ro(^i) and 
ro(^2) for system (Pq), with xo < ^1 < ^2- Let's call the intersection points among r 
and such level lines Pi _ , Pi^+ , with reference to £1, and P2- , P2,+ , with reference to 
£2, with 

P2,- -< Pi,- ^ P^ Pl,+ ^ ^2,+ . 

Analogously, when we consider two orbits r^(/ii) and r^(/i2) for system (P^), with 
Xfj. < hi < /i2, we name the intersection points among r and these level lines 
Qi- )Qi,+7 with reference to hi, and Q2- ,Q2,+, with reference to /i2, with 

Q2,- -< Qi,- ^ Q ^ Qi,+ -< Q2,+ ■ 
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Then the two annuh Ap and Aq turn out to be hnked together if 

P2,- -< Pi,- ^ Q2,- ^ Qi,- ^ Pi,+ -< P2,+ ^ Qi,+ -< Q2,+ . 

Proof of Theorem M.B. Consistently with the notation introduced above, we indicate 
the level lines filling Ap by ro(^), with £ e [^1,^2], for some Xo < ^1 < ^2, so that 

^p= U ToW- 

Analogously, we'll denote by V^{h), with h G [hi,h2\, for some x^i < ^1 < ^2, the 
level lines filling Aq, so that we can write 

^Q= U r^W- 

h\<.h<h2 

By construction, such annular regions turn out to be invariant for the dynamical 
systems generated by {Eq) and (-E^), respectively. We consider now the two regions 
in which each annulus is cut by the line r (passing through P and Q) and we call 
such sets ^p, ^p, A*Q and A^q, in order to have Ap = Ap U Ap and Aq = Aq U Aq, 
where the sets denoted by t are the "upper" ones and the sets denoted by b are the 
"lower" ones, with respect to the line r. We'll also name TZ'^ the rectangular region in 
which ^p and Aq meet and analogously we'll denote by 7^* the rectangular region 
belonging to the intersection between ^p and AU (see Figure 4). 




Figure 4. For the two linked annuli in the picture (one around P and the other 
around Q), we have drawn in different colors the upper and lower parts (with respect 
to the dashed line r), as well as the intersection regions 7?.* and TZ'^ between them. 
As a guideline for the proof, we recall that Ap is the annulus around P, having as 
its inner and outer boundaries the energy level lines ro(^i) and Tq{£2). Similarly, 
Aq is the annulus around Q, having as inner and outer boundaries the energy level 
lines Tf^{hi) and r^(/i2). 
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Let 

mi > 2 and m2 > 1 

be two fixed integers (the case mi = 1 and m2 > 2 can be treated in a similar manner 
and therefore is omitted). 

As a first step, we are interested in the solutions of system (£"0) starting from 
and crossing A^p at least mi times. After having performed the rototraslation of the 
plane that brings the origin to the point P and makes the x-axis coincide with 
the line r, whose equations are 




= {x — 2) cosu + — f ) sina; 

sinu; + — f ) cosa;, 



where a; := arctan(|), it is possible to use the Priifer transformation and intro- 
duce generalized polar coordinates, so that we can express the solution C(")^) — 
{x{-, z),y{-, z)) of system (£'0) with initial point in 2; = ixo,yo) £ «4.p through the 

radial coordinate p{t, z) and the angular coordinate 6{t, z). Therefore we can assume 
that ^^(0, z) G [— TT, 0]. For any t G [0, ro] and z G Ap, let's introduce also the rotation 
number, that is, the quantity 

roto(M)- 



27r 

that indicates the normalized angular displacement along the orbit of system (Eq) 
starting at z, during the time-interval [0,t]. The continuous dependence of the solu- 
tions from the initial data implies that the function (t, z) ^ 6{t, z) is continuous, as 
well as the map {t, z) 1— roto(t, z). From the definition of rotation number and by the 
star-shapedness of the level lines of Sq with respect to the point P, we have that for 
every z G To{£) the following properties hold: 

Vj G Z : roto(t, z) = j ^t = jto(^), 

Vj G Z : J < roto(t, z) < j + I ^ j To{i) < t < {j + 1) To{i) 

(if the annuli were not star-shaped, the inference " " would still be true). Al- 
though we have implicitly assumed that ii < i < £2 , such properties hold for every 
£>Xo- 

Observe that, thanks to the fact that the time-map Tq is strictly increasing, we know 
that To(-^i) < To (€2)- We shall use this condition to show that a twist property for the 
rotation number holds for sufficiently large time-intervals. Indeed, we claim that if 
we choose a switching time ro > a, where 

_ (mi + 3+|)To(4)ro(f2) 
ro{i2) - roiii) ' 

then, for any path 7 : [0, 1] Ap, with 7(0) G ro(-^i) and 7(1) G ro(-^2), the following 
interval inclusion is fulfilled: 



[9{ro, 7(1)), 9{ro, 7(0))] D [27rn*, 2n{n* + mi) - n], for some n* = n*{ro) G N. (2.2) 
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To check our claim, at first we note that for a path 7(5) as above, it holds that 
roto(t,7(0)) > [t/To{^i)\ and roto(t,7(l)) < \t/To{l2)'\, for every t > and so 

roto(t, 7(0)) - roto(t, 7(1)) > t """^t^T > 0. 

To (^1) To (£2) 

Hence, for t > a, with a defined as in (12.11) . we obtain 

roto(t, 7(0)) > mi + 1 + i + roto(t, 7(1)) , 

which, in turns, implies 

e{t, 7(0)) - e{t, 7(1)) > 27r(mi + 1), Vt > a. 

Therefore, recalling the bound 27r([t/ro(£2)l - f) < ^(^,7(1)) < 27r [t/ro(£2)l , the 
interval inclusion fl2.2p is achieved, for 

To 



n* = ri*(ro) : = 



(2.3) 



This proves our claim. 

By the continuity of the composite mapping [0,1] 3 s ^— roto(ro, 7(5)), it follows 
now that 

{e{ro, 7(s)), s G [0, 1]} 3 [27in*, 27r(n* + mi - 1) + tt]. 

As a consequence, by the Bolzano theorem, there exist mi pairwise disjoint maximal 
intervals [t/, t/'] C [0, 1], for i = 0, . . . , mi — 1, such that 

6'(ro,7(s)) G [2?™* + 2TTi, 2Tm* + tt + 27ri], Vs G i = 0, . . . , mi - 1, 

and 0{rQ, 'jiti')) = 2™* + 27ri, as well as 9{rQ, 'yit/')) = 27m* + vr + 27ri. Setting now 

7^l := and 7^2 := 7^^ 

we orientate such rectangular regions by choosing 

T^lieft ■= ^1 n ro(^i) and 7^-„^,, := 7^l n ro(£2), 

as well as 

'^2,ieft ■= ^2 n r^(/ii) and n^^right ■■= ^2 n r^(/i2) 

(see the caption of Figure 4, as a reminder for the corresponding sets). 
Introducing at last the mi nonempty and pairwise disjoint compact sets 

Hi := {z e Ap-. 6{ro, z) G [27rn* + 2Txi, 27m* + tt + 27ri]}, i = 0, . . . , mi - 1, 

we are ready to prove that 

{Hi, (j)o) :7Zi^7Z2, 2 = 0,...,mi-l, (2.4) 

where we recall that 0o is the Poincare map associated to system {Eq). Indeed, 
let's take a path 7 : [0,1] — > TZi, with 7(0) G T^ijeft 7(1) G Tli^right- For 
To > a and fixing i G {0,...,mi — 1}, we know that there exists a sub-interval 
[W, U"] C [0, 1], such that 7(t) G and 0o(7(t)) G ^p, Vt G [t/, t/']. Noting now that 
r^(0o(7(^/))) < ^1 and r^(0o(7(^/'))) — ^2, it follows that there exists a sub-interval 



16 M. PIREDDU, F. ZANOLIN 

[t*,tT] C [U'X] such that 0o(7(t)) e 7^2, Vt G and 0o(7(t:)) e ^2;^/*, as 

well as 00 (7 (^r)) ^ '^2, right- Therefore condition (12.41) is fulfilled. 

Let's turn to system {E^). This time we focus our attention on the solutions of 
such system starting from and crossing J\!q at least 1112 times. Similarly as before, 
we assume to have performed a rototraslation of the plane that makes the x-axis 
coincide with the line r and that brings the origin to the point Q, so that we can 
express the solution C(')^) of system {E^) with starting point in w G through 
polar coordinates (p, 6*). In particular it holds that 6'(0,w) G [0,7r]. For any t G 
[0, r^] = [0, T — ro] and w G ^g, the rotation number is now defined as 

~e{t,w)-~e{Q,w) 

ioi^[t,w) := . 

/TT 

Since the time-map is strictly increasing, it follows that T^{hi) < T^(/i2)- We claim 
that choosing a swithching time > /?, with 

o _ (^2 + 3 + |)r^(/zi)r^(/i2) 

then, for any path a : [0,1] Aq, with cr(0) G r^(/;,i) and cr(l) G r^(/i2), the 
following interval inclusion is satisfied: 

[e{r^, (t(1)), 9{r^„ a(0))] 3 [7r(2r2** + 1), 27r(n** + 7712)], for some n** = n**{r^) G N. 

(2.6) 

The claim can be proved with arguments analogous to the previous ones and therefore 
its verification is omitted. The nonnegative integer n** has now to be chosen as 

-jh] ■ (2.7) 

By (12.61) and the continuity of the composite mapping [0, 1] 9 s 1-^ rot^(r^, (t(s)), 
it follows that 

{e{r^„ cr{s)), se[0,l]}D [2™** + vr, 27r(n** + m^)]. 

As a consequence, Bolzano theorem ensures the existence of m2 pairwise disjoint 
maximal intervals [s/, s/'] C [0, 1], for i = 0, . . . , 777.2 — 1, such that 

6'(r^, a{s)) G [27m** + vr + 27ri, 27m** + 27c + 27ri], Vs G [s/, s/'], i = 0, . . . , 7773 - 1, 

and 6'(r^, o-(s/)) = 27777** + tt + 27r7, as well as 9{r^, cr(s/')) = 27r77** + 27r + 2777. 

For TZi and 7^2 as above and introducing the 7712 nonempty, compact and pairwise 
disjoint sets 

/Cj := {w G Aq : 6'(r^, 77;) G [2?™** + tt + 27rz, 2?™** + 27r + 27rz]}, z = 0, . . . , 7772 - 1, 

we are in position to check that 

(/C„ 0^) : n2^nu 7 = 0, . . . , m2 - 1, (2.8) 

where (p^ is the Poincare map associated to system (-E^). Indeed, taking a path 
(T : [0, 1] — » 7^.2, with a{0) G T^2,ieft and (t(1) G T^2, rights (3 and for any 7 G 
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{0, . . . , m2 — 1} fixed, there exists a sub-interval [s/, s/'] C [0, 1], such that cr(t) G /Cj 
and G Vt G [s/, s/']. Since ro(0^(a(s/))) < £i and ro(0^(a(s/'))) > ^2, 

there exists a sub-interval [s*,s**] C [s/, s/'] such that (j)^,{(j{t)) G T^i, Vt G [s*,s**] 
and 0^(cr(s*)) G T^ijeft^ ^^^^ 4>ii{(^{s**)) G T^i right- Therefore condition fl2.8p is 
proved. 

The stretching properties in (12.41) and (12.81) allow to apply Lemma 13.11 and the 
thesis follows immediately. □ 

We observe that in our proof we have chosen as TZi the "lower" set TZ'' and as 7^2 the 
"upper" set 7^*. However, since the orbits of both systems (Eq) and (-E^) are closed, 
the same argument also works (by slightly modifying some constants, if needed) if we 
choose TZi = TZ* and 7^2 = 7?.^. 



3. Topological tools and remarks about chaotic dynamics 

In this last section, we briefiy recall the topological tools that we have employed 
along the paper. They are taken, with minor variants, from |32j and are based on 
some previous works [30l [311 [33l [35]. We start with some terminology: 
Given a metric space Z, we define a path in Z as a continuous map 7 : M D [0, 1] ^ Z 
(instead of [0, 1] one could take any compact interval [sq, Si]). A sub-path cr of 7 is the 
restriction of 7 to a compact sub-interval of its domain. An arc is the homeomorphic 
image of the compact interval [0, 1]. By a generalized rectangle we mean a set 7?. C Z 
which is homeomorphic to the unit square Q := [0, 1]^ C M^. 

If 7^ is a generalized rectangle and h : Q ^ h{Q) =TZis a homeomorphism defining 
it, we call contour dTZ of TZ the set 

m ■= h{dQ), 

where dQ is the usual boundary of the unit square. Notice that i!}7l is well defined, as 
it is independent of the choice of the homeomorphism h. We name oriented rectangle 
the pair 

n := (7^,7^"), 

where 7?. C Z is a generalized rectangle and 

7^-:=7^^,^,u7^-^,, 

is the union of two disjoint compact arcs TZ^^j^^ , 7^7ig/it — ''^'^^ ^^^^ '^^^^ ^^ft 
the right sides of 7Z~ . Once that 7Z~ is fixed, we can also define TZ'^ as the closure of 
{}7Z \ IZ^ . In particular, we set 

7?+ — 7?+ I ) 7?+ 

where TZ'^own ^"^^ '^up t"^o disjoint arcs. In the usual applications, the ambient 
space Z is just the Euclidean plane and the generalized rectangles are bounded by 
some orbit-segments, possibly associated to different systems. 

The central concept in our approach is that of "stretching along the paths" : 
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Definition 3.1. Suppose that ip : Z D ^ Z is a. map defined on a set and 
let X := {X,X~) and Y := {Y,Y^) be oriented rectangles of a metric space Z. Let 
/C C X n be a compact set. We say tliat (/C, ip) stretches X to Y along the paths 
and write 

if the following conditions hold: 

• ip is continuous on /C ; 

• for every path 7 : [0, 1] — X such that 7(0) G Xj^j-^ and 7(1) G -^^gM (o'^ 
7(0) G X~gi.^^ and 7(1) G Xj^j^), there exists a sub-interval C [0, 1] such 
that 

7(t)G/C, v^(7(t))Gr, wte[t',t"] 

and, moreover, 'ip{'-){t')) and ip{l{t")) belong to different components of Y~ . 
In the special case in which /C = X, we simply write 

^ : X^F. 

We are now in position to introduce the result about the existence of chaotic dy- 
namics for linked twist maps that we apply in the proof of Theorem ll.2[ 

Lemma 3.1. Let Z he a metric space, let ^ : Z D Dq, Z and \1' : Z D Dx^ ^ Z he 
continuous maps and letTZi := {TZi,TZi~), 7I2 := (7^2, "^2") he oriented rectangles in 
Z. Suppose that the following conditions are satisfied: 

there exist mi > 1 pairwise disjoint compact sets Hi , . . . ,Ti.mi ^ "^i ^ 
such that (Hi, $) : T^i^^ 7^2, for z = 1, . . . , mi ; 

(-ff^) there exist m2 > 1 pairwise disjoint compact sets /Ci , . . . , /Cma ^ "7^2 H 

such that (/Cj, : TZ2^^ T^i, for i = 1, . . . ,m2 ■ 

If at least one hetween mi and m2 is greater or equal than 2, then the composite map 
'?/':= \i/ o $ induces chaotic dynamics on mi x m2 symbols in the set 

H* := [j where H[j ■.= H,n<^-\lCj). 

i— 1, .. .,m-^ 
j = l,...,m2 

Moreover, for each sequence of mi x m2 symhols 

s = (s„)„ = {pn, qn)n G {1, . . . , mi}^ X {1, . . . , ma}^, 
there exists a compact connected set Cg C H'p^^ g^ with 

Cs n 7^+,,^„ ^0, n 7^+„p ^ 

and such that, for every w E Cg , there exists a sequence {yn)n with y^ = w and 

yn e H'p^^q^ , i){yn) = Vn+i , Vn > 0. 
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Note that in Theorem 11.21 we haven't used the last part of the previous lemma. 
Therefore it would be possible to improve our main result on Lotka-Volterra systems 
by adding information about the existence of continua of points for which we have a 
control on the forward itineraries. 

We end the paper by trying to clarify the relationship between the concept of 
chaos expressed in Definition 11.11 and other ones available in the literature, with 
special reference to the semiconjugation with the Bernoulli shift and the density 
of the periodic points. Before giving the next lemmas we just recall some basic 
definitions. We denote by = {1, . . . , m}^ the set of the two-sided sequences of m 
symbols. Analogously, by S+ = {1, . . . , m}^ we mean the set of one-sided sequences 
of m symbols (where N is the set of nonnegative integers). Introducing the standard 
distance 

d{s', s") := J2 , where s' = {s',U , = (^D^a , (3-1) 

for I = Z and for I = N, we make and compact metric spaces. 

Our first result is quite classical. Indeed, the same conclusions of Lemma 13.21 have 
been obtained by several different authors dealing with similar situations (see, for 
instance, ^Ql Theorem 3]). Nonetheless, we give a detailed proof for the reader's 
convenience. We also thank Duccio Papini for useful discussions about this topic. 

Lemma 3.2. Let {Z, dz) he a metric space, ip : Z D Z be a map which induces 

chaotic dynamics on m > 2 symbols in a set V C , relatively to (/Ci, . . . , /Cm) 
(according to Definition \l.l\) . Assume also that ip is continuous and injective on 

m 

i=i 

Define the nonempty compact set 

J:= fl r^(/C). (3.2) 

j=-oo 

Then I is invariant for ip andip\x is semiconjugate to the two-sided m-shift, through 
a continuous surjection g : I ^ as in (11.61) . Moreover, the counterimage through 
g of any k-periodic sequence in contains at least one k-periodic point of I. 

Proof. First of all, we observe that w G X if and only if there exists a full orbit, that 
is, a two-sided itinerary, {wi)i^z such that Wq = w and ip{wi^i) = Wi E JC for every 
i G Z. By the assumptions on ip coming from Definition 11.11 and standard properties 
of compact sets, it follows immediately that X is a nonempty compact set such that 
^(X)=X. 

As a next step, we introduce a function gi , associating to any G X its corresponding 
full orbit 
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that is, the sequence of points of the set X defined by 

Wi := ip^{w), Wi G Z, 

with the usual convention ip'^ = Idjc and = ip. The injectivity of ip imphes that 
the map 

is well defined. Recalling that the )C/s are pairwise disjoint, we know that for every 
term Wi of there exists a unique label 

Si = Si{wi), with Si G {1, . . . ,m}, 

such that Wi G /C^- . Hence, the map 

92-Si,t-^ (si)igz e T^m (3.3) 

is also well defined. Thus, if we set 

by Definition ll.il we obtain a surjective map that makes the diagram ( 11. 6p commute. 
Moreover, the inverse image through g of any A;-periodic sequence in contains at 
least one /c-periodic point of I. 

Finally, we show that g is continuous. To such aim we call 

e := min dz{ICi,)Cj) > 

and note that the map ip\x : I = X is a homeomorphism. Therefore the 

following property holds: 

Vra G N, 3 5 = 5„ > : for u,v el, with s„ := (Mi)igz , := {vi)iez , 

dz{u,v) < 6 =^ dz{ui,Vi) <e, W\i\ <n =^ Si{ui) = Si{vi), W\i\ < n. 

From this fact and the choice of the distance in (see (13. ip ). the continuity of g is 
easily proved. □ 

The next result is just a more precise version of Lemma [3. 2 [ with reference to the 
periodic points. 

Lemma 3.3. Under the same assumptions of Lemma \3.2l there exists a compact 
invariant set A C X such that is semiconjugate to the two-sided m-shift, through 
the continuous surjection g\x (where g : I ^ is like in (11. 6p ). The set of the 
periodic points of tplj is dense in A and, moreover, the counterimage through g of any 
k-periodic sequence in contains at least one k-periodic point of A. 

Proof. For X defined as in (13. 2p . we consider the subset V made of the periodic points 
of ip\i , that is, 

V := {w e 1 : 3k > 1, ^\w) = w} 

and define 

A := V. 

Since ipiV) = V, it follows that tp{A) = A. From the last statement in Lemma [3.21 we 
also find that g(V) coincides with the subset of made by the two-sided periodic 
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sequences of m symbols, which is dense in • This latter fact implies the surjectivity 
of : A — > . The remaining properties are a straightforward consequence of the 
corresponding ones in Lemma 13.21 □ 

In both the above results, we have required the injectivity of the map ip. This is 
not a heavy restriction in the applications to ODEs (like those in Section [1] and in 
j), where ip is the Poincare map. Anyway, in the abstract setting of the metric 



spaces and with reference to Definition 11.11 it is still possible to obtain some suitable 
versions of Lemma 13.21 and Lemma 13.31 for a map ip which is continuous on /C, but 
without the assumption of injectivity. Such task can be accomplished in different 
ways. A first possibility is that of considering a semiconjugation with the Bernoulli 
shift on m symbols for one-sided sequences. Indeed, any initial point w uniquely 
determines the forward sequence := (^/'*(w))jgN , to which we can associate a well 
determined sequence of symbols (si)jgN £ S+ as in the definition of g2 in (13. 3p . 
Another possibility is that of considering, in place of g, only the map g2, that is, to 
associate to any full orbit {wi)i^z, with 'ip{wi) = Wj+i , the sequence (sj)jgz ^ 
such that Wi G /C^. . This second approach is followed by Lani-Wayda and Srzednicki 
in [22] , where the authors obtain a variant of Lemma 13.31 for a set A replaced by the 
closure (denoted by T) of the set of the terms of the periodic sequences in /C which 
are full orbits of (and project through g2 onto a periodic sequence of S^). Hence, 
in this case, the commutative diagram (11.61) becomes 



92 



92 



where V'^((wj)iGz) := (wii+i)is 



4. Conclusion 

It is a pleasure and a honor to have the possibility to dedicate our work to the 
memory of Professor Andrzej Lasota, who gave fundamental contributions both in 
the area of chaotic dynamics and in the study of mathematical models for biological 
systems. 
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